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Abstract 

We study three different (/-analogues of the harmonic numbers. As applications, 
we present some generating functions involving number theoretical functions and 
give the ^-generalization of Gosper's exponential generating function of harmonic 
numbers. We involve also the g-gamma and g-digamma function. 
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1 Introduction 

The harmonic numbers are defined as 

n i 

H n = J2T, H :=0. (1) 

k=i K 

Our aim is to find some g-analogues of these numbers. We start from some 
basic identities satisfied by the harmonic numbers. For example, 
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H n = -s{n + 1,2), (2) 
n! 

where s(n + 1, 2) is a Stirling number of the first kind [6]. 
Next we introduce the most basic notions of g-calculus. Let 

and 

[n] fl ! = [nUn-l] a ...[l] g . 

Define also (x; q) n = (l-x)(l-qx) ■ ■ ■ {l-q n ~ l x) and (x; q)^ = lim n ^ 00 (x; q) n . 
Then the (/-binomial coefficient with parameter n and k is 



r n\ [n) q \ (q;q), 



(n>k> 0). 



\kj q [k] q \[n-k] q \ (q;q) k (q;q) n - k 
According to ([I])-© we define three class of g-harmonic numbers as 



H n,q 



H; 



n,q 



k=i We 

r s,(n + l,2), 



n 



fe=i \ k / u 



-1 



[*] s 



i7. 



n,q 



u- 1 



(4) 
(5) 
(6) 
(7) 



(Let Hq :— for i — 1, 2, 3, 4.) Here s g (n + 1, 2) is a g-Stirling number of the 
first kind PQ. 

The first definition appears in [4"fTT] , for example. We point out that the first 
and second definitions are the same but the third and the fourth are different 
from each other - although all of them tend to H n as q — > 1. That differs 
from the others is obvious because of the presence of hi(g). 

In what follows we deduce some identities and give applications involving these 
numbers. 
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2 Identities involving and H% q 

2. 1 Number theoretical results 



In this section we point out that the H\ g-harmonic numbers are connected to 
the divisor function. Moreover, the products of the Riemann zeta function and 
polylogarithms are Dirichlet generating functions of some interesting number 
theoretical "polynomials" . 

By definitions, 

Cauchy's product gives that 

E ff »/ = r^EA ( 8 ) 



n>l n>\ 



We need the notion of Lambert series [7J p. 257]. In general, a Lambert series 
has the form 

n 

n>i 1 y 

where a n is any suitable sequence. This connection implies some interesting 
results. We cite a useful theorem of [7J Theorem 307]: if 

n>l n>\ 



then 



n 

Hi) = E = E M r ' 

n>l 1 y n>l 



holds if and only if 

C(s)/(s) = 

where 



n>l 



is the Riemann zeta function. 
Let us apply this and (JSJ): 



E W = E = 7^-^(9) (9) 
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with A simple transformation shows that 



n>l \d\n J 

Let us introduce the notation 

A(x,n) = £x d . 

d\n 

In special, 

A(l,n) = d(n) = a (n), 
where d(n) is the divisor function and 

<7x(n) = £<f. 

d\n 

It is worth to realize the symmetry between the power and the base in A(x, n) 
and a x (n). 

Substituting x = 1 in (jUJ) we have 
Proposition 1 

J2H n>q q n =J2l n d(n). 

n>l n>l 

So the generating function of the divisor function involves the (/-harmonic 
numbers. More generally, with the help of (jUJ) and f JTOj) and the definition of 
A(x, n) one can write 

E H n , q (qx) n = £ q n A(x,n). 

n>l 1 Q X n>l 

Again, by © 

n 

n>l 1 y n>l 

Hence we choose and b n = A(x,n) and apply the theorem cited above: 



Li s (x) = E^7> and 9{s) = Y, Ai ' '" 



„>i n „>i n 



then 

Proposition 2 



C^Li^HE^'" 



»>i ra 
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Here Li s (x) is the well-known polylogarithm function. In special Li s (l) = £(s) 
so we get the known sum 

d{n) 



C\s) = £ 



n>l nS 



2.2 A recursion for H\ 



Since 



we may think that 



k=l 



q 

holds for some sequence a*-. This is really true but, sadly, does not have a 
simple form. See the table below. 

We shall need the notion of the g-Seidel matrix [2]. Given a sequence a n , the 
g-Seidel matrix is associated to the double sequence o\ given by the recurrence 



a„ = a n (n> 0) 

k „n„k—l i „fc— 1 



< = «" i + <;l (n>0,*>l). 

In addition, a° is called the initial sequence and Oq the final sequence of the 
g-Seidel matrix. Then the identity 

t (?) * (12) 

k=0 



connects the initial and the final sequence. 
Define the generating functions of a° and 



and 



a(x) = a °n xn i a ( x ) — a o xn i 

n>0 n>0 



a(x) = a °rri ' A ( x ) = »<> 



n>0 l n \l- n>0 y n \l m 



A proposition given in |2] states that these functions are related by the fol- 
lowing equations: 



5 



*(*) = £<fc^— . ( 13 ) 

n > l X > <?Jn+l 

A(x) = e g (x)A(x), (14) 
where 



>0 W 



is the g-analogue of the exponential function [5]. We introduce the notation 
Egf(a n ) and Gf(a n ) for the exponential and ordinary generating function of 
a n , respectively. 

To reach our aim posed in ( 1111) . our approach is as follows. Let the final se- 
quence b n = H*. We determine the initial sequence a° = a n . Then Egf(6 n ) = 
Egf(Hng) = e g Egf(a n ). And, to get Egf(a n ) we determine a n by using ([8]) and 

Gf(M ^ Gf«) = £ ^- = £ ^7^—- (15) 

1 ~~ X n>l 1 _ 9 „>1 I s ) Q)n+1 

From this equation a n can be determined. (Note that a = 0.) 
Proposition 3 We have 



n 



H n,q = S ( ) G * 



fc = l V'"/ <j 

where the sequence is determined recursively by 



£i V fc -V g H 1-9 



Proof. The denominator of the right hand side of ( fl5l) is 

1 11 1 (q n qx; q) c 



(x;q) n+1 l-x(qx;q) n 1-x (qx;q) c 
The g-binomial theorem [51 Section 1.3] states that 



(az;g)oo ^ (a; k 



Applying this to f fT6|) . 



1 (q n qx;q) 00 _ 1 ^ (g";g)fc ^ fc 



1-x (qx-^q)^ l-xf^ (q;q) k 



(16) 
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Thus (TT5l) becomes 



i-g)E = E «nx n I E ^P^M* 



Let 



R (g n ; g)fc fc 
(?;?)* 



for short. Then 



for all n and A;. Moreover, 



-Bfc,n = g 



'n + k - 1 N 



(1 - g) E -&-n = E E B k ,nX k . (17) 

n>l 1 g n>0 \fc>0 / 

If we write the sums term by term, we get 
«o(-Bo,o + B li0 x + B 2 fiX 2 H ) + aiar^-Bo,! + + -B2,i^ 2 H ) + •■■ = 

a B 0>0 + x(a B lfi + OiS ,i) + x 2 (a B 2fi + a 1 B ljl + a 2 B 0i2 ) H 

Comparing the coefficients here and the left hand side of ( jTTl) . we have 

n 1 - g 

E dkB n - ktk = • 

/c=o 1 y 

Note that - bacause of (|T5l) - ao must be zero. Remember also that is the 
initial sequence of our g-Seidel matrix, so ( 1121) gives 

H L = E Q «*■ (18) 

This is our proposition. 

Remark. It is worth to present the first terms of the sequence a n : 



a 
en 



a 3 



a 5 



= 0, 
1, 



a 2 = -- 



q 2 + q - 1 



g+ 1 



g 5 + g 4 



g + 1 



g 2 + g + 1 



04 = — - 



9 i 8 

q + g 



2g 5 + q 2 + q - 1 



g 3 + q 2 + q + 1 
g 14 + g 13 — g 10 — g 9 — g 8 



g 7 + g 6 + g 5 — g 2 — g 



g 4 + g 3 + g 2 + g + 1 
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_q 20 + g 19 - g 16 - 2g 14 + q 12 + q 11 + q w + q 9 - 2q 7 - q 5 + q 2 + q - 1 
6 q 5 + q 4 + q 3 + q 2 + q + 1 

It seemed to be interesting to give a simple formula for the nominator. How- 
ever, one can easily see that 

a k -)• - — f as q -> 1 (fc = 1,2,3,4,5,6). 



According to ([3]), this is plausible for all ak- 

As a consequence of ( 1T41) and ( IT8l) . we have the next connection: 

Egf« g ) = e,Egf(a n ). 

This is a curious version of Gosper's identity (}2"Tj) involving the exponential 
generating function of the harmonics. 



2.3 The case H 2 



To close the case of H\ , we remark also that 



H L = H l, ( 19 ) 



(see (j5J) for the definitions). The (/-Stirling numbers of the first kind are defined 
recursively by [TJ p. 155] 

s q (n + 1, k) = s q (n, k - 1) + [n] q s q (n, k), (20) 

and s q (0, 0) = 1, s q (n, 0) = when n > 0. 
These relations imply 

H ^ = ^ Sq{n + 1,2) = W Sq{n ' l) + W^W^ 2) ' 



hence 



Hl„ = Hi, + — — s q (n, 1) = Hl_ x 



n-l] q \ 



Wig! 



Then (EDI) gives (EE 
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3 Identities involving H% 

3.1 A q-analogue of Gosper's result 



The exponential generating function of the harmonic numbers is deduced by 
Gosper 



n 



n>l 



X 

J2 H n —r = xe x 2 F 2 



1 


1 1 










V 


2 2 





(21) 



Here 



ai, a 2 , 



E 



„>0 ( & l)n(&2)n ■ ■ • (&s)n ™' 

is the hypergeometric function with parameters a» and 6-,-, and 

(a) k = a(a + l)(a + 2) ■ • • (a + k - 1) 
under the agreement (a)o = 1. 

The (/-version of the hypergeometric function is called basic hypergeometric 
function and defined as [5] 



E 



[ ai, ei2, 

\ 6l, &2, ••• 

(ai; g)n(a2; (a r ; 9)7 



-1)V 



(g;?)t 



(bi;q) n (b 2 ;q)n---(b s ;q)n 
One may see from (122]) that 

(& + l)th term x(-q h ) l+s ~ T (1 - a 1 q k ) •••(! — a r q k 



(22) 



fcth term 



1 - q k+1 (1 - hq k ) ••■(!- 6 s g A 



Now we derive the (/-analogue of Gosper's result fl2~T|) with respecto to 
Equations (EJ) and (TH1) give that 



Egf(fl? .) = e g Egf g 



Hence it is enough to determine the sum 



E? 

k>l 



-l) k+1 X k 
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We would like to express this function as a basic hypergeometric series. Thus 
k should run from 0, so we write 



t„ q [k+l], [fc+l],!' 



whence 



(k + l)th term _ (—xq)q k (1 — gg fc ) 2 



Since s 
state the following 



kth term ~ LI \ — q k+1 (1 — q 2 q k ) 2 

r = 2 in our case, the place of (—1)™ is indifferent. Finally we can 



Proposition 4 VKe have 



EH 3 
n,q 

n>l 



X 



n 



1 - q)xe q (x) 2 (f>2 
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q; -qxj 









3.2 The Hockey Stick Theorem - a q-analogue 



It is natural to ask, what is the recursion satisfied by H 3 ? If H 3 = H^_ ± + 



i,q ' n,q n — l,q 

would be true, we knew that ifj? „ = H}. „. But this relation does not hold. 

[n\ q ' n,q n,q 

In fact, the next recursion is valid. 



Proposition 5 For H 3 



3 , „ l n ~ ^iq 



We see that the limit of q — —, , is 1 — — -. So the standard recursion 
for harmonic numbers holds just asymptotically. To prove the proposition, 
we need the following statement. The standard (not q-) version can be found 
at the webpage |http: / /binomial.csueastbay.edu/IdentitiesNamed.htm[ Cata- 
log #: 3400001. 

Proposition 6 (Hockey Stick Theorem) 



n—j 

EC- 1 ) 

fc=l 




The name comes from the fact that the summands and the sum has a shape 
in the Pascal triangle like a hockey stick. 
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We prove the (/-analogue version. 

Proposition 7 (Hockey Stick Theorem — a g-analogue) 



'n-Y 



Proof. Write the sum term by term: 

( n 



+ 



n 



1)" Y ■ ' I " I + (-l) n - j+l q(*) 
\ n - 1 J q 

The binomial coefficients are rewritten with the recursion [5] 

n — l\ In — l\ 



n 
. k , 



k 



k - 1 



(j+i) 



n — 1 



+ 



q 



n-1 

v 3 



q\ 2 



,j+3 1 



7 9 
71 — 1 



+ 



71—1 



J+2 



9 

/ 72 — 



+ 



n — 1 



j + 3 \j + 2 J 

\ J / q \ J / q 



J + 2 \j + 1/ 



+ 



-1)^3 



In — 


1 + 


In — 


3 


\n — 


' q 


\n — 





+ 



.1)^+^(2) 







'n-r 

,72-1. 



Realize that the first members in the square brackets with the (/-power coeffi- 
cients are cancelled by the second member in the next square bracket: 



/j+k+i\ I n — 1 \ 

,r "\j+k)-* ' V*). 



0. 



This is true for all k — 1, 2, . . . ,n — j — 1. 

Now we are ready to prove Proposition [51 A consequence of ([6]): 

/ n — l\ \ g( 2 ) 



^n.o n n-l, 



n 
,1, 



1 lq) [1], 



+ 



n \ in — 1 



[2]« 



'n-r 

, , 



(1-9) 



'n-l N 
. 1 , 



'n-Y 
v 3 , 



q 



(2) 



[3], 



+ • ■ 



In — 1 



(2) 



[1]< 



+ 



1 v^Twjp], 
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The first sum is exactly the left hand side of the g-Hockey Stick Theorem with 
j — 1, up to a minus sign. The second sum equals — 1 [TJ p. 153]. Thus 



H ^- H "-^ =1 -w q { n 1 1 ) q +q ~ 1 - 

So we have the desired result. 



4 Identity involving 



The digamma function is defined as 

r'(:r) 



where T(x) is the Euler gamma function. Latter satisfies the functional equa- 
tion 

T(x + 1) = xT(x). (23) 

The harmonic numbers are connected to the digamma function (as one imme- 
diately sees by the logarithmic derivative of ( I2"3"j) ): 

H n = iP(n + 1) + 7- 

Here 7 = —"0(1) is the Euler-Mascheroni constant. Our goal is to find the 
g-analogue of this formula. 

Let us start from the definition of the g-gamma function (see [5]): 

r « (l) = fefe (1 -^- 
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The g-digamma function is simply the logarithmic derivative of T q (x) as in 
the ordinary case ( Pl8][T0] ): 



ip q (x) = - ln(l -q) + ln(g) ^ - 



q 



n+x 

n>0 ' 



It is known that T q (x) satisfies the next q- version of (j23 

T q (x + 1) = [x] q Tg(x). 
Then taking logarithm and then the derivative, we get 

ip q (x + l) = -^ln([x] 9 ) + if> q {x). 

Since 



dx 9 q x — 1 ' 
the previous equation is rewritten recursively as 

^ — 1 q x 1 — 1 q — 1 

So definition ([7]) seems to be correct. 

On the other hand, we may define the g-analogue of the Euler-Mascheroni 
constant as 

n 

lq = -^(1) = ln(l - q) - ln(g) £ , 

n>l 1 y 

or more simply, 
With these we have 
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